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Abstract. The self-consistent Theory of Finite Fermi Systems based on the Energy Density Functional by 
Fayans et al. with the set DF3-a of parameters fixed previously is used to calculate three kinds of quadrupole 
moments. At first, we examined systematically quadrupole moments of odd neighbors of semi-magic lead and 
tin isotopes and N = 50, N = 82 isotones. Second, we found quadrupole moments of the first 2 + states in the 
same two chains of isotopes. Finally, we evaluated quadrupole moments of odd-odd nuclei neighboring to double 
magic ones. Reasonable agreement with available experimental data has been obtained. Predictions are made for 
quadrupole moments of nuclei in the vicinity of unstable magic nuclei 
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1 Introduction 

As the result of quick development of experimental tech- 
niques in nuclear physics, the bulk data on nuclear static 
moments has become very extensive and comprehensive 
Q], thus creating a challenge to nuclear theory. First of all, 
it concerns the nuclei distant from the /J-decay stability val- 
ley which are often close to the drip lines and are of great 
interest to nuclear astrophysics. For this reason, a theoreti- 
cal approach used for describing such nuclei should have a 
high predictive power. The self-consistent Theory of Finite 
Fermi Systems (TFFS) |2| based on the EDF by Fayans et 
al. is one of such approaches. 

A good description of the quadrupole II4I5H and magnetic 1 6 7 
moments of odd semi-magic nuclei has been achieved within 
this approach. For quadrupole moments, we use a new DF3- 
a version [8 | of the original DF3 functional [3] which was 
employed in calculations of magnetic moments. It differs 
from the DF3 version only in the spin-orbit parameters 
H, x' and the effective tensor force. The DF3-a functional 
is characterized by a rather strong effective tensor force. 

In these calculations, the so-called "single-quasiparticle 
approximation" has been used, where one quasiparticle in 
the fixed state A = (n, I, j, m) with the energy sa is added 
to the even-even core. According to the TFFS Q, a quasi- 
particle differs from a particle of the single-particle model 
in two respects. First, it possesses the local charge e q and, 
second, the core is polarized due to the interaction between 
the particle and the core nucleons via the Landau-Migdal 
(LM) amplitude. In other words, the quasiparticle possesses 
the effective charge e e ff caused by the polarizability of the 
core, which is found by solving the TFFS equations. In 
the many-particle Shell Model, a similar quantity is intro- 
duced as a phenomenological parameter which describes 
polarizability of the core consisting of outside nucleons. It 
should be noted that for this series of problems within the 
scope of the mean-field theory, the self-consistent TFFS is 
similar to the HF-QRPA approach. 



Recently, quadrupole moments of the first 2 + state in 
even lead and tin isotopes have been found IflOl . again 
within the self-consistent TFFS. For this problem which is 
evidently beyond the mean-field theory, the TFFS results 
are significantly different from the QRPA ones. 

In this paper, we review briefly the results of the cited 
references on quadrupole moments of odd nuclei and of 
the 2 + states in even-even ones and add some new calcula- 
tions for odd unstable nuclei. In addition, we include here 
some results of flU for the quadrupole moments of odd- 
odd nuclei neighboring to the double magic ones. 



2 Brief calculation scheme 



The calculation scheme of the self-consistent TFFS based 
on the EDF method by Fayans et al. is described in de- 
tail in Ref . |4j . Here we write down only several formulas 
which are necessary for understanding main ingredients of 
the approach. The EDF method by Fayans et al. (3) is a 
generalization for superfluid finite systems of the original 
Kohn-Sham EDF method |[T2"1 . In this method, the ground 
state energy of a nucleus is considered as a functional of 
normal and anomalous densities, 



J £[p„(r), p p (r), v„(r), Vpir)]^ 



(1) 



Within the TFFS, the static quadrupole moment Qx of 
an odd nucleus with the odd nucleon in the state A can be 
found in terms of the diagonal matrix element (A\V(u> 
! A) of the effective field V in the static external field 
Vo = V167r/5r 2 i20- In systems with pairing correlations, 
equation for the effective field can be written in a compact 
form as 
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V(co) = e q V (co) + TA{aS)V{aS), 



(2) 
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where all the terms are matrices. In the standard TFFS no- 
tation [9], we have: 



(3) 



Table 1. Quadrupole moments Q (e b) of odd-neutron nuclei in 
the state A. 
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(5) 



where £, Mi, and so on stand for integrals over e of the 
products of different combinations of the Green function 
G(e) and two Gor'kov functios F^'(e) and F (2) (e). They 
can be found in J9]- 

Isotopic indices in Eqs. d3]|5]l are omitted for brevity. 
The explicit form of the above equations is written down 
for the case of the electric (f-even) symmetry we deal with. 
In Eq. (|4), T is the usual LM amplitude, 



T = 



6 2 8 
dp 2 ' 



(6) 



is the density-dependent effective pairing interaction, 

6 2 & 



5v 2 



(7) 



and the amplitudes T ^ = stand for the mixed second 
derivatives, 

= — . (8) 

8p8v 

In the case of volume pairing, one has J ro ^ = 0, whereas 
for the case of surface pairing we deal the amplitude is 
non-zero and should be taken into account when Eqs. (|3]|5]l 
are solved. As the analysis of Ref. [4| shows, the compo- 
nent V of the vector V, as a rule, dominates. However, the 
fields d\,d 2 also contribute, and sometimes significantly, 
to the value of Qx- 

In this article the TFFS equations are solved in the 
self-consistent basis obtained within the EDF method with 
the functional DF3-a. Thus, the same set of parameters 
has been used to calculate the single particle scheme and, 
according to Eqs. (161718b . the effective interactions in the 
TFFS equations. We consider the surface kind of pairing 
as motivated by our previous research [4|, see also ab ini- 
tio arguments in Ref. [ 1 3 1 . 



3 Quadrupole moments of odd semi- and 
near-magic nuclei 



nucleus 


X 




<2,h 


6Q 


J9 Ca 


1*4/2 


0.036(7) 


+0.040 


0.004 






0.040(6) 




0.000 


41 Ca 


l./*7/2 


-0.090(2) 


-0.078 


0.012 






-0.066(2) 




-0.012 






-0.080(8) 




0.002 


85 Kr 


139/2 


+0.443(3) 


+0.507 


0.064 


87 Kr 


2d 5/2 


-0.30(3) 


-0.355 


-0.06 


87 Sr 


109/2 


+0.33(2) 


+0.335 


0.01 


89 Sr 


2d 5/2 


-0.271(9) 


-0.245 


-0.026 


89 Zr 


1^9/2 


+0.28(10) 


+0.262 


-0.02 


9I Zr 


2^5/2 


-0.176(3) 


-0.195 


-0.019 






(-)0.257(13) 




0.062 






-0.206(10) 




0.011 


10< >Sn 


2^5/2 


+0.31(10) 


+0.250 


-0.06 


m Sn 




+0.18(9) 


+0.029 


-0.13 


115 Sn 




0.26(3) 


+0.377 


0.12 


ll9 Sn 


2d V2 


0.094(11) 


-0.035 


-0.129 






-0.065(5) 




0.030 






-0.061(3) 




0.026 


121 Sn 


2d 3/2 


-0.02(2) 


+0.063 


0.08 


135 Xe 


2d 3/2 


+0.214(7) 


+0.217 


0.003 


137 Xe 


2/7/2 


-0.48(2) 


-0.376 


0.10 


137 Ba 


2d 3 /2 


+0.245(4) 


+0.254 


0.009 


139 Ba 


2/7/2 


-0.573(13) 


-0.445 


0.128 


141 Nd 


2d i/2 


+0.32(13) 


+0.289 


-0.03 


143 Nd 


2/7/2 


-0.61(2) 


-0.518 


0.09 


143 Sm 


2^3/2 


+0.4(2) 


+ 0.296 


0.1 


145 Sm 


2/7/2 


-0.60(7) 


-0.537 


0.06 


197p b 


3P3/2 


-0.08(17) 


+0.195 


0.27 




3P3/2 


+0.08(9) 


+0.272 


0.19 


201 nu 

Pb 


2/5/2 


-0.01(4) 


+0.137 


0.15 


203 ra. 

Pb 


2/5/2 


+0.10(5) 


+0.284 


0.18 


205 m. 

Pb 


2/5/2 


+0.23(4) 


+0.336 


0.09 


209 m, 

Pb 


2^9/2 


-0.3(2) 


-0.264 


0. 1 


211m. 

Pb 


2^9/2 


+U.U9(6) 


-O.ZSi 


n in 
-U.3 / 


U3 Sn 


1/l l*l/2 


0.41(4) 


-0.776 


-0.37 




0.48(5) 




-0.30 


113 Sn 


1/! l*l/2 


0.38(6) 


-0.703 


-0.32 


117 o 

U/ Sn 


lh U/2 


-0.42(5) 


-0.593 


-0.17 


U9 Sn 


lh' n/2 


0.21(2) 


-0.469 


-0.25 


121 Sn 


1/j* 
1 "ll/2 


-0.14(3) 


-0.293 


-0.15 


l23 Sn 


1^11/2 


+0.03(4) 


-0.123 


-0.15 


l25 Sn 


1^11/2 


+0.1(2) 


+0.039 


-0.1 


135 Xe 


lfl U/2 


+0.62(2) 


+0.504 


0.12 


I37 Ba 


lh' n/2 


+0.78(9) 


+0.588 


-0.19 


,47 Gd 


U l3/2 


-0.73(7) 


-0.791 


-0.06 


191p b 


U 1V2 


+0.085(5) 


+0.0004 


-0.085 


193p b 


1! 13/2 


+0.195(10) 


+0.335 


0.140 


195p b 


1! 13/2 


+0.306(15) 


+0.689 


0.383 


197p b 


l! 'r 3/ 2 


+0.38(2) 


+0.980 


0.60 


205p b 


h 'l*3/2 


0.30(5) 


+0.665 


0.37 



The final expression for the quadrupole moment of an odd 
nucleus is as follows 1191 141 : 

er = (4-fi)v p A (9) 

where u\, v\ are the Bogolyubov coefficients and 

Va = f V(r)R 2 nlj (r)r 2 dr. (10) 



For odd neighbors of a magic nucleus the "Bogolyubov" 
factor in (0 reduces to 1 for a particle state and to -1 for 
a hole one, see also lfl5l . If the odd nucleon belongs to 
the superfluid component, the factor (u 2 - v 2 ) in Eq. (0 
becomes non-trivial. It changes permanently depending on 
the state A and the nucleus under consideration. This factor 
determines the sign of the quadrupole moment. It depends 
essentially on values of the single-particle basis energies 
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sa reckoned from the chemical potential yu as we have 

(u 2 A -v 2 A ) = (e A -n)/E A , (11) 



Ea = yje 2 + d 2 . Keeping in mind such sensitivity, we found 
this quantity for a given odd nucleus (Z, + 1) or (Z+ 1, N), 
N, Z even, with taking into account the blocking effect in 
the pairing problem [14| putting the odd nucleon to the 
state A under consideration. For the Va value in Eq. (01 we 
used the half-sum of these values in two neighboring even 
nuclei. 

The results of the calculations are presented in Tables 1 
and 2 which contain odd-neutron and odd-proton nuclei re- 
spectively with known experimental quadrupole moments 
(our predictions for odd nuclei with unknown quadrupole 
moments see in [4 5 1). One can see that the theoretical sign 
of the quadrupole moment is correct in all cases when the 
sign of the experimental moment is known. This permits 
to use our predictions to determine the sign when it is un- 
known. Several rather strong disagreements with the ex- 
perimental data for high-j levels l/zn/2 in Sn isotopes and 
h'13/2 in Pb isotopes originate from their too distant posi- 
tions from the Fermi level, see [5 1, where it was found that 
the Q values depend strongly on the single-particle level 
structure. It follows mainly from Eq. (11). 

To evaluate the agreement with experiment quantita- 
tively, we calculate the mean theory-experiment difference 

V^=^Z(e; h -er p ) 2 , d2) 

with obvious notation. On average, the agreement, can be 
considered as reasonable. For 42 quadrupole moments of 
odd-neutron nuclei, the average disagreement between the- 
ory and experiment is not so small, 

j(8Q)L =0.189 eb. 
However, it is concentrated mainly in 15 intruder states for 

which we have -\j (c50 2 ms [intruder] = 0.269 e b. For the 
rest of 27 "normal" states, the disagreement is rather mod- 
erate yj (6Q) 2 ms [normal] = 0.125 e b. For protons, agree- 
ment is worse. The rms deviation is >/(*0L =0.254 eb. 
The main contribution to this deviation comes from In and 
Sb isotopes, odd neighbors of even tin nuclei. It is the result 
of too strong quadrupole field V lup (r) for the DF3-a func- 
tional \4\. For neutrons, this drawback is partially hidden 
with multiplying by the Bogolyubov factor, but for protons 
it appears to the full extent. For more detailed discussion, 
see 0. 

For odd-neutron neighbors of even N=50 isotones, the 
proton-subsystem is superfluid and the neutron Bogolyubov 
factor in Eq. (0 is ±1. In this case, agreement with the 

data is almost perfect, 

V^rms =0.041 e b. The situa- 
tion is similar for odd-neutron neighbors of even isotones 
with N=82. Again, agreement with the data is rather good, 

V(£0L=O.O93 eb. 

For the major part of nuclei in Tables 1 and 2, neighbor- 
ing to double-magic ones, let us call them "near-magic", 
the quality of agreement is rather good. Therefore, we hope 
to predict reasonably the quadrupole moment values for 
such nuclei including strongly proton- or neutron-rich ones. 
These predictions are presented in Table 3. 



Table 2. Quadrupole moments Q (b) of odd-proton nuclei in the 
state A. 



nucl. 


A 


Qexp 


2,h 


SQ 




ld 3/ 2 


0.0585(6) 


0.069 


0.010 


41 Sc 


I/7/2 


-0.156(3) 


-0.139 


0.017 






0.120(6) 




-0.019 






0.168(8) 




0.029 


87 Rb 


2^3/2 


+0.134(1) 


+0.132 


-0.002 






+0.138(1) 




-0.006 


105 In 


109/2 


+0.83(5) 


+0.833 


0.00 


107 In 


109/2 


+0.81(5) 


+0.976 


0.17 


109 In 


109/2 


+0.84(3) 


+1.113 


0.27 


m In 


109/2 


+0.80(2) 


+ 1.165 


0.36 


1I3 In 


109/2 


+0.80(4) 


+ 1.117 


0.32 


1I5 In 


109/2 


+0.81(5) 


+ 1.034 


0.22 






0.58(9) 




0.45 


117 In 


109/2 


+0.829(10) 


+0.965 


0.136 


119 In 


109/2 


+0.854(7) 


+0.909 


0.055 


121 In 


109/2 


+0.814(11) 


+0.833 


0.019 


123 In 


109/2 


+0.757(9) 


+0.743 


-0.014 


i 25ln 


109/2 


+0.71(4) 


+0.663 


-0.05 


127 In 


109/2 


+0.59(3) 


+0.550 


-0.04 


U5 Sb 


2d 5 /2 


-0.36(6) 


-0.882 


-0.52 


U9 Sb 


2d 5/2 


-0.37(6) 


-0.766 


-0.40 


121 Sb 


2d 5/ 2 


-0.36(4) 


-0.721 


-0.36 






-0.45(3) 




-0.27 


123 Sb 


107/2 


-0.49(5) 


-0.739 


-0.25 


137 Cs 


107/2 


+0.051(1) 


-0.031 


-0.080 


139 La 


107/2 


+0.20(1) 


+0.103 


-0.10 


141 Pr 


2d 5/2 


-0.077(6) 


-0.120 


-0.043 






-0.059(4) 




-0.061 


145 Eu 


2d 5 /2 


+0.29(2) 


+0.156 


-0.13 


205 j j 


3d* /2 


+0.74(15) 


+0.227 


-0.51 


203 Bi 


lh 9/2 


-0.93(7) 


-1.323 


-0.39 






-0.68(6) 




-0.64 


205 Bi 


lh<)/2 


-0.81(3) 


-0.945 


-0.14 






-0.59(4) 




-0.36 


207 Bi 


lh 9/2 


-0.76(2) 


-0.454 


0.31 






-0.55(4) 




0.10 






-0.60(11) 




0.15 


209 Bi 


lh<)/2 


-0.516(15) 


-0.342 


0.18 






-0.37(3) 




0.03 






-0.55(1) 




0.21 






-0.77(1) 




0.43 






-0.40(5) 




0.06 






-0.39(3) 




0.05 


213 Bi 


lh 9/2 


-0.83(5) 


-0.508 


0.32 






-0.60(5) 




0.09 



As it was mentioned in Introduction, the core polar- 
izability by the quadrupole external field is characterized 
directly by the effective quadrupole charges, which are de- 
fined naturally within TFFS as e'^ = V^ n /(V%) A |fl6l 
In Tables 1 and 2, there are only two nucleus, 209 Bi and 
2()9 Pb, with a double-magic core. In this case, the polariz- 
ability is relatively moderate, e^ ff = 1.4, e"^ = 0.9. In nu- 
clei with unfilled neutron shell, it becomes much stronger, 
e e ff — 3-^6 |01 ■ The reason is rather obvious. Indeed, for the 
case of positive parity field Vo, virtual transitions inside the 
unfilled shell begin to contribute in such nuclei and small 
energy denominators appear in the propagator X", Eq. (5), 
playing the main role in Eq. (fJJ for the problem under con- 
sideration. This enhances the neutron response to the field 
Vq and, via the strong LM neutron-proton interaction am- 
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Table 3. Predictions for quadrupole moments Q (e b) of odd 
near-magic nuclei. 





j 


T, 

1/2 


CI , 


n 


!Ni 2V 


7/2+ 


204.7 ms 


-0.26 


- 




3/2- 


35.6 h 


-0.17 


- 


2gNi 49 


(9/2) + 


128 ms 


0.20 


- 


SNiso 


(5/2 + ) 


635 ns 


-0.12 




[01 o_ 
50 bn 51 


(5/2) + 


1.7 s 


-0.21 




50 ^"81 


(3/2 + ) 


56 s 


0.10 


-0.04(8) 




7/2- 


1.46 s 


-0.17 




207pu 
g 2 "125 


(1/2)- 


stable 







27C028 


7/2" 


17.53 h 


0.31 




29CU28 


3/2" 


196.3 ms 


-0.20 




29CU50 


(3/2-) 


188 ms 


-0.13 




49 iri50 


(9/2) + 


3 s 


-0.35 




131 Tn 


(9/2) + 


0.28 s 


0.28 




133 cu 
51 oD 8 2 


(7/2) + 


2.34 m 


-0.23 




207 Tl 
81 il 126 


1/2+ 


4.77 m 









Fig. 1. g 2 order corrections to the mass operator in magic nuclei. 
The circles with one wavy line in the first term are the phonon 
creation amplitudes g. The second term is the phonon tadpole. 




Fig. 2. Matrix element M LL in the form of Eq. (16). 



phonon: 

M LL = V°(r)df>G(r,r,s)dr—, (13) 



plitude T np , the proton response as well. The results in Ta- 
ble 2 for the chain 203,205,209gj sn ow how the polarizability 
grows with increase of the number of neutron holes. Keep- 
ing in mind this physics, one can represent the effective 
charges as e^ ff = 1 + e p x , e" ff = where efjj is the pure 
polarizability charge. To separate contributions of the un- 
filled shells and core nucleons explicitly, one can divide 
the Hilbert space of the QRPA equations (f2]i to the "va- 
lent" and subsidiary ones and carry out the corresponding 
renormalization procedure ifTTl . 



4 Quadrupole moments of the first 2 + 
states in Sn and Pb isotopes 

Account for the phonon coupling (PC) is the direct way 
to generalize the standard nuclear theory. As a rule, the 
so-called g 2 approximation is used where g is the phonon 
creation amplitude. However, almost all of these general- 
izations did not take into account all the g 1 terms, limit- 
ing themselves with the mass operator pole diagrams only, 
see the first diagram in Fig. 1, where diagrams for the 
mass operator are displayed. The second diagram repre- 
sents the sum of all g 1 non-pole diagrams usually called 
the tadpole. The problem of consistent consideration of all 
g 1 terms including tadpoles was analyzed firstly in the ar- 
ticle by Khodel [ 1 8 1 . The method developed was applied 
to magic nuclei, mainly for ground state nuclear charac- 
teristics, within the self-consistent TFFS (2). It was found 
that, as a rule, the tadpole contributions in magic nuclei 
are noticeable and are often of opposite sign as compared 
with those of the pole terms. The first attempts to include 
phonon tadpole effects for nuclei with pairing were re- 
cently made in Refs. |[T9l, lEUl and IflOl. 

According to Ref. [18], in the g 2 approximation, the 
matrix element Mn for a static moment of the excited state 
(phonon) with the orbital angular moments L in a static 
external field V° is determined in terms of the change of 
the one-particle Green function (GF) in the field of this 



6%G = 6 L (Gg L G) = G(e)g L G(e + u L )g L G{e) (14) 
+G(e)g L G{e - u) L )g L G(e) + G{e)5 L g L G{e), 

where gi is the amplitude for the production of the L phonon 
with the energy u>l and 5igL is the variation of gi in the 
field of other L phonon. This quantity is the main part of 
the phonon tadpole in Fig. 1 . After some transformations 
of these expressions one can obtain 

M LL = V°Gg L Gg L G + V°A6 L g L . (15) 

It is convenient to transform this expression in such a 
way that the effective field V, Eq. ©, appears instead of 
the external field V°. After regrouping terms in Eq. (15[ 
and in the integral equation for &l9l, for details, see Refs. 
B2I20I10I . we obtain the ultimate expression, 

M LL = VGg L Gg L G + VAS L TAg L , (16) 

which is illustrated in Fig. 2. It contains now the effective 
field V, Eq. (0, instead of V° and the quantity 5lT in the 
second term which denotes the variation of the effective ph 
interaction T in the field of the L phonon. For the density 
dependent TFFS effective interaction Tip), the following 
ansatz can be readily obtained [ 1 8j2) : 

&F 

SlT{y) = — p£(r)F iM (n), (17) 
dp 

where p\ = Ag^ is the transition density for the L phonon 
excitation. The first term of Eq. (16) coincides with the 
result of Refs. [21 22] while the second one, with the 5lT 
quantity, is a generalization to take into account all the g 2 
terms. 

All the above equations can be readily modified for 
such processes as the transition between the excited states 
L and L' in the external field V°(u> = u>u - cjl) or the ex- 
citation of the two-phonon state L + L' in the external field 
V°(u> - (±>u + o>l). The static moment case corresponds to 

LO — 0, ix)u — Ll>l- 

This approach for magic nuclei has been generalized 
for non-magic ones in 11231 1 01 . Then eight matrix elements 
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Table 4. Quadrupole moments Q (e b) of the first 2 + states in Sn 
and Pb isotopes. 



Fig. 3. Matrix elements for M L [ and M,£ for non-magic nuclei 




100 104 108 112 116 120 124 128 132 
A 

Fig. 4. Quadrupole moments of the first 2+ excited states in even 
Sn isotopes. 



instead of one in Fig. 2 should be considered, two of them 
are shown in Fig. 3. 

It is necessary to compare this expression with the QRPA 
approach. Here we mean the usual scheme [24 1 which uses 
the QRPA wave functions for the matrix element between 
two excited states. In Ref. [24] the expression forthe B(E2) 
quantity has been derived using the bare external field V° 
and the QRPA wave functions without the pp and hh-channels. 
The analytical expression for the sum of the eight above- 
mentioned matrix elements consists of two parts. The first 
part coincides with the corresponding formula in 0241141 
with one important correction, which is the first generaliza- 
tion of the QRPA approach. Namely, instead of the external 
field y°, which does not depend on the frequency, the ef- 
fective field V appears, which depends in general on the 
frequency a> = ll>i + u' L . The second part of the sum is new 
and describes the contribution of the ground state corre- 
lations (GSC), the so-called backward-going diagrams, to 
the first diagrams of Fig. 3 with the integrals of three GF's 
("triangle"). This is the second generalization. We calcu- 
late the contribution of such correlations separately. The 
terms with 5lT and 5lT^, Fig. 3, are the third generaliza- 
tion of the QRPA approach. Note that these terms are also 
absent in Refs. H22I2 1 1261271251 . The main difference of 
our approach from the calculations in Refs. 122 25 2 61271 
is the self-consistency on the (Q)RPA level and absence of 
any phenomenological or fitted parameters. 

We calculated the quadrupole moments of the first 2 + 
states in the tin and lead isotopes in the /l-representation 
with self-consistent single-particle wave functions <px ob- 
tained within the EDF method of Ref. [3] with the func- 
tional DF3-a (HI. A spherical box of the radius R- 16 fm 
is used to simulate the single-particle continuum. We ex- 
amined the dependence of the results on the cut-off energy 
£max and have found that the value of £' max =100 MeV en- 
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— 


-0.03 


-0.001 


132 Sn 


0.04 


— 


0.05 


0.015 


lj4 Sn 


-0.01 




0.00 


-0.001 


I90p b 


-0.92 




-0.30 


-0.008 


192p b 


-1.15 




-0.38 


-0.008 


194p b 


-1.31 




-0.44 


-0.008 


l%p b 


-1.26 




-0.42 


-0.008 


198 Pb 


-1.05 




-0.35 


-0.008 


200p b 
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0.09 
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0.06 


-0.002 
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0.05 


-0.7(3) 
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0.043 



sures 1% accuracy. To calculate the quantities V and gi, 
the results of Ref. H have been used where all the calcu- 
lations were performed in the coordinate representation us- 
ing the same self-consistent DF3-a basis as in the present 
calculation of the matrix element Mn. Thus, the single- 
particle continuum is taken into account adequately in the 
present calculations. The contribution to the Q values of 
the term with 5T in Fig. 3, Eq. (1171 1. turned out to be rather 
small, (-(0.01 +0.03) e b). However, there are cases where 
these corrections are comparable with the total Q(2^) value 
when the proton and neutron values almost compensate 
each other, e.g. in I18 Sn and 122 Sn nuclei. The term with 
5lT^ contains the anomalous analogs of the correspond- 
ing quantities in Eq. (11711 , 

The results are given in Table 4 and Fig. 4. Except for 
1I2 Sn and 208 Pb nuclei, we obtained a reasonable agree- 
ment with experimental data [1 ]. The contribution of the 
GSC term turned out to be large. Often it is more than 
50+60% of all triangle contributions (column Q(GSC=0)). 
The usual QRPA (GSC=0 and V = V°), see the last column 
in Table 4, results in considerably less Q values. 
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5 Quadrupole moments of odd-odd near 
magic nuclei 

As the odd-odd nuclei are more complicated objects than 
the odd ones, we consider here only the near-magic odd- 
odd nuclei. Within the above-described self-consistent EDF 
approach, we calculated the ground state quadrupole mo- 
ments of odd-odd near-magic nuclei with the use of the ap- 
proximation disregarding the interaction between two odd 
particles. This simple approximation can be checked in a 
pure phenomenological way [28 11] and it turned out rea- 
sonable. Within this approximation, the problem is reduced 
to calculations of quadrupole moments of corresponding 
odd nuclei. 

Indeed, if we neglect the interaction between two quasi- 
particles, the quadrupole moment of the odd-odd nucleus 
with the spin / is as follows: 

Qi =< // | V + V" | // >, (18) 

where Wu = £(pi<f2 < j\ m \h m i I H >> for the particle- 
particle case. Here ip\ is the single-particle wave function 
with the quantum numbers 1 = A\ = (n\, ji,h,m{). Then 
the expression for ground state quadrupole moment of the 
odd-odd near-magic nucleus has the form: 



Table 5. 

nuclei. 



Quadrupole moments Q (e b) of odd-odd near-magic 



Q, = (2I - l)|5o_ 7 7 |( 1)' 



jp I jn 
I jp 2 



jn I jp 
I jn 2 



+j„+I+2 



(19) 



where 3j-symbol cj = 2 j(2 j - 1) 1/2 [(2 j + 3)(2 j + 2) x 
x(2j + l)2j]~ , Q p and Q" are the quadrupole moments 
of corresponding odd nuclei. Similar formulae can be eas- 
ily obtained for the hole-hole, particle-hole and hole-particle 
cases. The details can be found in the Poster article by 
Voitenkov et al. ifTTI at the given conference. 

The results of the calculations are presented in Table 
5. In column Q e ff we show the results for the odd-odd nu- 
clei under consideration obtained with the effective charges 
e P e ^ = 2, e"yy = 1 in order to compare them with this well- 
known phenomenological description. We see that there 
are only three experimental Q values for the long-living 
nuclei and our approach describes them rather satisfactory. 
Other nuclei are short-living ones and our results give reli- 
able predictions for their Q values. 



6 Conclusion 

Quadrupole moments of odd neighbors of semi-magic lead 
and tin isotopes and N = 50, /V = 82 isotones are calcu- 
lated within the self-consistent TFFS based on the Energy 
Density Functional by Fayans et al. with the DF3-a param- 
eters fixed previously. The same approach has been used 
to calculate quadrupole moments of the first 2 + state in tin 
and lead isotopes as well as the moments of near-magic 
odd-odd nuclei. 

For the quadrupole moments of odd and odd-odd near 
magic nuclei a good agreement with the experiment has 
been obtained. For the case of semi-magic nuclei a reason- 
able agreement with experiment for the quadrupole mo- 
ments has been obtained for the most part of nuclei con- 
sidered. In this case when the odd particle belongs to the 
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superfluid subsystem, the Bogolyubov factor (u A - v\) = 
(sa —fi)/EA comes to the quadrupole moment value, in ad- 
dition to the matrix element of the effective field Va- This 
factor makes the quadrupole moment value very sensitive 
to calculation accuracy of the single-particle energy ea of 
the state under consideration, especially near the Fermi 
surface as the quantity Qa vanishes at e\ = /U. For such 
a situation, influence of the coupling of single-particle de- 
grees of freedom with phonons, see [6 20 5 1, should be es- 
pecially important. 

For the quadrupole moments of the first 2 + states, we 
have obtained a noticeable difference from the traditional 
QRPA approach. In particular, new terms with 6lF and 
&lF^ appear, which contain the density derivatives of both 
the ph and pp effective interactions. In the problem un- 
der consideration, their contribution turned out to be rather 
small, as a rule. However, for consistency, these terms should 
be included. Except for the 112 Sn and 208 Pb cases, a rea- 
sonable agreement has been obtained with the experiment 
available. Using the self-consistent method which contains 
no newly adjusted parameters we have also predicted the 
values of quadrupole moments of the first 2 + states in sev- 
eral unstable lead and tin isotopes including the 100 Sn and 
132 Sn nuclei. An unexpectedly large contribution of ground 
state correlations to the Q(2\) values is found. A non-trivial 
dependence of the quadrupole moments of the first 2 + states 
on the neutron excess is found which can be traced to the 
negative proton contributions. A similar behavior could prob- 
ably be present in other isotope chains. 

Quadrupole moments of unstable nuclei including those 
near the exotic 100 Sn and 132 Sn 56 78 Ni are also predicted, 
which should be of special interest. 
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